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All scalar master integrals (Mis) for massive 2-loop QED Bhabha scattering are 
identified. The 2- and 3-point Mis have been calculated in terms of harmonic 
polylogarithms with the differential equation method. The calculation of 4-point 
Mis is underway. We sketch some alternative methods which help to solve (mainly) 
singularities of some Mis. 



1 Introduction 

Bhabha scattering, e + e~ — » e + e~(nj), is used in many accelerators for the 
determination of luminosity. The NLO corrections in the electroweak model 
are known 1 . The lack of knowledge of the complete NNLO corrections in 
massive QED has been one of the major sources of the theoretical error at 
LEP 2 , the other two are light fermion initial state pair production and the 
hadronic vacuum polarization. Several calculations have been already done to 
get massive NNLO QED results for Bhabha scattering 3 ' 4,5,6 . The first step in 
calculations is to identify Mis. For genuine vertices the task has been completed 
in 3 . The additional 3-point Mis coming from boxes and, more involved, the 
4-point Mis themselves are presented for the first time at this conference in 
the first part of the talk. The second major step includes the evaluation of 
Mis. This process is underway. Two powerful methods are in use. One is 
based on a bottom-up approach using differential equations 8 . The other uses 
Mellin-Barnes representations for the calculation of each MI separately. The 
(semi)-analytical results are already known for some of the more complicated 
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cases 5 . Further, there is a general numerical method to calculate Mis at fixed 
kinematical Euclidean points 9 . 

For our complete set of Mis and the current status of calculations, see 6 
and the web page 7 . Here we focus on a presentation of the second part of the 
talk, namely some alternative cross-checks of analytical results. 

2 Some cross-checks of analytical results 

2.1 Algebraic relations between IR-divergent Mis 

The method of determining Mis, which is realized in the package IdSolver 10 , is 
based on the Laporta-Remiddi (LR) algorithm 11 : it determines and solves an 
appropriate set of algebraic equations with integration by parts 12 and Lorentz 
invariance 13 identities. The result is a file, sometimes huge in size, which 
involves relations among Mis with different powers of propagators and irre- 
ducible numerators. Such relations can be used to fix singularities of purely 
IR divergent Mis. An example is shown in Fig. 1. 




Figure 1: Relation between a vertex OBJ (here V512m2 in terminology of 6,7 ), the same 
vertex with a dotted line OBJd (here V512m2d), and with irreducible numerator (fcipi). 

The diagrams OBJ and OBJd are UV finite, but IR divergent. The corre- 
sponding diagram with an irreducible numerator is finite. If we expand the 
known Mis (second row in Fig. 1) into series in e = (4 — d)/2 with known 
coefficients Fi/e n \ and make an ansatz for the unknowns, 

OBJd(s) = S_ 2 (s)i + £_i(s)~ + B (s) + ■ ■ ■ , (1) 

and similarly for OBJ, we get from Fig. 1 relations among the coefficients. 
Here, e.g. (a, 6, c, d are known singularities of the appropriate simpler Mis): 

0=4 [aF 3 (e) + bF 4 (e) + cF 5 (e) + dF 6 (e) + F 2 (e)£?_ 2 ( S )] , (2) 
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which allows to determine B_ 2 (s) = l/(4s), in agreement with 6 . The crucial 
point is that there is an additional factor 1/e in front of OBJd in comparison 
to OBJ. To get e.g. £?_i(s) or A- 2 (s) (analogous coefficient of OBJ), another 
equation with a different, independent numerator would have to be used in 
addition. 

2.2 Exact subloop integration 

The diagram in Fig. 2 has a massless UV divergent subloop. 




d d k 1 d d k 2 

[k 2 2 ][(ki + k 2 -Pi) 2 ][kl +m»][(fci + p2f 



Figure 2: A diagram with a massless propagator subloop. 



This allows to perform the subloop integration (over k 2 ) analytically, resulting 
in a Co-like function where one of the denominators appears with power e. 
Using a Feynman parameter representation, after integration over k\ and one 
of the Feynman parameters, we get (t = (pi + p 2 ) 2 ): 

d r(l-e) 2 r(e) T(2e) 
16r(l + e) 2 r(2-2e) 1 + j r(2 + e) ™ m - 

This generates a singularity in I num at a; = 0, which can be treated by the 
following subtraction (F(x,y) = [(1 — x)(l — y) 2 + xy^]~ 2e ): 

dxdy x- 1+e (l ~ X y- 2e 
o ([(1 ~ - y) 2 m 2 + xyt} 2 ) 2 ' 

i r r(e)r(2-2 £ ) i 

T(2 - e) 1 - 4e + reg 



(m 2 ) 2c 

7 reg - / dxx- l+e (l-.,)> Jl / ^//[Flr./yl ■ r(().//l] 



l /.l 
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The remaining integrations in I reg can be performed analytically or numerically 
after the e-expansion: 

I reg = dx(l ~ x y** e -*Mi-.) f d y [ln/( "' V) ln/( °' y)] 
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E(- 2e ) 7 



n=0 



J^ln^- 1 f(x,y)la k f(x,0) 



k=0 



he + I 2 e 2 



In this way both the singularities and also regular terms can be obtained. We 
have checked that they coincide with our results in 6 and 7 for the MI V411m2. 

2.3 Subtracting a counter term 

The diagram V411ml, drawn in bold lines in Fig. 3, is UV divergent. 




Villml 



(!) -kO ^ x> 

Iji + ll -\\ +i[i + 2 + 2lDg(l- a; )-bg( a; )] 



Figure 3: The divergences of the diagram V411ml can be obtained by renormalization theory 
methods. 



Let us introduce a subtraction as in Fig. 3. This subtraction is actually a 
counter term for the subdivergence in the MS scheme. It is known 14 that 
after such a renormalization of subdivergences the divergences of the result are 
polynomials in dimensionful parameters. Since the diagram is dimensionless, 
they are just constants and we can replace massless lines by massive ones 
and set the external momentum to zero, in this way avoiding spurious IR 
divergences. This replacement is shown in the second row of Fig. 3 and allows 
to calculate the singular part of the diagram. For this, we also have to add the 
subtraction at arbitrary s, x = ^fijl"^ ! (second square bracket in Fig. 3). 
The final result agrees with 6 and 7 for the MI V411ml. 



2.4 Expressing a subtracted subdiagram by a dispersion relation 

The diagram V413md contains an IR singular subloop when the dotted line 
becomes on-shell. After integrating out the UV divergent 2-point subloop, we 
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get 



Ivmmd = I jj— 2 _ — ^r^B [{k + q) 2 ,m 2 ,m 2 ]. (3) 



By subtracting and adding the 2-point subloop as shown schematically in Fig. 
4 we get a finite vertex type integral, plus a product of 2-point functions. 




~ k + q~0~ % O 

Figure 4: The IR divergent 2-loop diagram V413md with a subtraction. 

The vertex integral may be further rewritten using a dispersion relation rep- 
resentation for the difference Bo[(k + q) 2 ,m 2 , m 2 ] — Bo(q 2 , m 2 ,m 2 ) under the 
integral. After all these preparations, the e expansion of the diagram may be 
determined. 

Methods like the ones presented in this section are also applicable for 4- 
point functions. 
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